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Iterative data reconstruction 



5 The present invention relates to the field of iterative data reconstruction, 

e.g. in computed tomography. In particular, the present invention relates to a method of 
performing an iterative data reconstruction, to an image processing device and to a 
computer program for an image processing device. 

10 

Iterative methods can be used for data reconstruction in different fields 
of data processing, such as nuclear science and computed tomography. In particular, 
iterative algorithms are applied for image reconstruction. The reconstruction process 
consists of alternating projection and backprojection steps. In order to efficiently use 

1 5 dedicated image reconstruction hardware, multiple projections/backprojections have to 
be carried out simultaneously. In known methods, for example, as described in 
IP 0 502 187 Bl or US 6,574,299 Bl, such image reconstruction still requires long 
processing times. In particular, the application of multiple projections/backprojections, 
which may be carried out simultaneously, may lead to a slow convergence, due to 

20 correlations between the simultaneously processed data and, consequently, to a low 
image quality. 

It is an object of the present invention to provide for an improved 
25 iterative data reconstruction. 

According to an exemplary embodiment of the present invention as set 
forth in claim 1, the above object may be solved by a method of performing an iterative 
data reconstruction, wherein projection data is determined from estimated data for a 
plurality of projections. Then, a difference between the estimated data and measured 
30 data is determined. Then, according to an aspect of the present invention, a filtering of 
the difference is performed, resulting in a filtered difference. Then, a backprojection is 
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performed by updating the estimated data by using the filtered difference. These steps 
may be performed iteratively. 

Due to the filtering step, the difference is adapted or manipulated before 
performing the backprojection. This may allow for an improved processing speed, i.e. 
5 for a reduced processing time. Furthermore, in case the method is applied, for example, 
for the reconstruction of images, an improved image quality may be achieved. 

According to another exemplary embodiment of the present invention as 
set forth in claim 2, the filtering is performed such that a mutual influence or reciprocal 
interaction of the plurality of projections is at least partially filtered out. In other words, 
10 according to this exemplary embodiment of the present invention, the difference is 
modified such that an influence on a projection caused by other projections is 
compensated for before backprojection. 

Advantageously, due to the introduction of such a filtering step, the 
method according to this exemplary embodiment of the present invention may be 
15 implemented efficiently, for example, on dedicated image reconstruction hardware and 
may allow to perform several projections/backprojections simultaneously, allowing for 
a reduced processing time. 

According to another exemplary embodiment of the present invention as 
set forth in claim 3, the method is based on the algebraic reconstruction technique 
20 (ART). 

Claims 5 to 7 provide for further exemplary embodiments of the method 
according to the present invention. 

According to another exemplary embodiment of the present invention as 
set forth in claim 8, an image processing device is provided, performing an iterative 
25 data reconstruction, for example, similar to the ART, including a filtering before a 
backprojection, allowing for a reduced processing time, while still allowing for a high 
reconstruction quality. 

The present invention also relates to a computer program, for example, 
for an image processing device, for performing an iterative data reconstruction, 
30 including a filtering step. The computer program according to the present invention is 
defined in claim 9. The computer program according to the present invention is 
preferably loaded into a working memory of a data processor. The data processor is 
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thus equipped to carry out the method of the invention. The computer program may be 
stored on a computer readable medium, such as a CD-ROM. The computer program 
may also be presented over a network, such as the Worldwide Web, and may be 
downloaded into the working memory of the data processor from such a network. 
5 It may be seen as the gist of an exemplary embodiment of the present 

invention that a filtering step is introduced into an iterative data reconstruction, such as 
ART or simultaneous ART (SART). ART is, for example, described in R. Gordon et al 
"Algebraic reconstruction techniques (ART) for three-dimensional electron microscopy 
and x-ray photography" J. Theor Biol. Vol. 29, pages 471 to 481, 1970, which is hereby 

10 incorporated by reference. SART is, for example, described in R.H. Andersen et al 9 

"Simultaneous algebraic reconstruction technique (SART)" Ultrasonic imaging, Vol. 6, 
pages 81 to 94, 1994, which is hereby incorporated by reference. The filtering 
performed according to the present invention allows to filter out influences caused by 
other projections onto the current projection, which allows to improve the quality of the 

15 data reconstruction, i.e. in case images are reconstructed, it allows for an improved 
image quality. Furthermore, according to the present invention, several projections / 
backprojections may be performed simultaneously, allowing for a high processing 
speed. 

These and other aspects of the present invention will become apparent 
20 from and elucidated with reference to the embodiments described hereinafter. 

Exemplary embodiments of the present invention will be described in the 
following, with reference to the following drawings: 



25 Fig. 1 shows a schematic representation of an image processing device 

according to an exemplary embodiment of the present invention, 
adapted to execute a method according to an exemplary 
embodiment of the present invention. 
Fig. 2 shows a thorax phantom, reconstructed with one iteration of ART 
30 withX=l. 

Fig. 3 shows a comparison of SART (left) and a reconstruction 

performed in accordance with an exemplary embodiment of the 
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present invention (right) for M = 32 (1 iteration, X = 1). Where 
M is the number of simultaneously processed views. 
Fig. 4 shows another comparison of S ART (left) and a reconstruction 
performed in accordance with an exemplary embodiment of the 
5 present invention (right) for M = 64 (1 iteration, X = 1). 

Fig. 1 depicts an exemplary embodiment of an image processing device 
according to the present invention, for executing an exemplary embodiment of a 
10 method in accordance with the present invention. The image processing device depicted 
in Fig. 1 comprises a central processing unit (CPU) or image processor 1 connected to a 
memory 2 for storing projection data and the data generated during the iterative data 
reconstruction. The image processor 1 may be connected to a plurality of input/output- 
network - or diagnosis - devices, such as an MR device, or a CT device. The image 
15 processor 1 is furthermore connected to a display device 4 (for example, to a computer 
monitor) for displaying information or images computed or adapted in the image 
processor 1. An operator may interact with the image processor 1 via a keyboard 5 
and/or other output devices which are not depicted in Fig. 1. 

The image processing device depicted in Fig. 1 may be operated on the 
20 basis of the ART reconstruction technique. 

The basic idea of ART used according to the present invention is to use a 
discrete representation / of the continuous object function and to calculate projection 
data from it. lis modified, if there if a difference between the calculated projection data 
and the measured data. 

Let the measured projection data p consist of X views pi*..,px. 
One iteration step &M> k + 1 consists of two operations: 
1: For a given view n(k) projection data p' is calculated from an 
estimated image I k and is compared with the measured data p n(k) . (projection) 

P ~ ^n{k)^k 

P n(k) denotes the projection operator for view n(k) . 
2: The estimated image is updated depending on the observed difference 
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between the measured and the calculated projection, resulting in a new estimate I M . 
(backprojection) 

B n(k) is the backprojection operator for view n(k) . 

5 n is the sequence in which the projection data from different views is 

processed (i.e. n : iVl— > {l,...,X}). X is a weighting factor that controls how much of the 
observed difference is backprojected into the new image. 

A drawback of ART is that the computational effort is fairly high. This 
may be overcome by using a dedicated image reconstruction hardware, such as a 
10 dedicated CPU which can calculate several projections/backprojections simultaneously. 
As one iteration step in ART consists of one projection/backprojection pair, ART has to 
be modified to use several projections simultaneously. This leads to the simultaneous 
algebraic reconstruction technique (SART), which may be used to operate the image 
processing device depicted in Fig. 1. 
15 In SART Mprojections/backprojections are processed simultaneously in 

each iteration step k h-> k + M : 

1 : Projection data p* . is calculated from an estimated image / k and 

compared with the measured data p n(k+J) for all j e [0 V ..,M -l]. (projection) 
p'j = P n{k+j) h V,e[0,...,M-l] 
20 Let A j bs X n(k+J y(p n(k+J) ~p\j ) 

2: The estimated image is updated depending on the observed difference 
between the measured and the calculated projection, resulting in a new estimate I k+M . 
(backproj ection) 



J M-l 



25 The factor 1/M in the backprojection step is due to the fact that 

projections from different angles partly contain the same information about the object. 
For example, all projections contain the dc value (overall average) of the object. While 
the factor 1/M is adequate for the dc component, it is too high for higher frequency 
components. This leads to slow convergence. 
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According to a preferred embodiment of the present invention, the image 
processing device depicted in Fig. 1 is adapted to perform the following method of 
operation, which is referred to in the following as filtered SART. 

According to this exemplary embodiment of the present invention, a 
5 filtering step is introduced, which is performed before backprojection. Advantageously, 
the filtering may be adapted such that a mutual influence of the plurality of projections 
may be filtered out, thus allowing for an improved image quality. In other words, the 
slow convergence due to correlations between the simultaneously presented data may 
be significantly improved and thus the image quality may be improved. According to an 
10 exemplary embodiment of the present invention, the filtered SART may be described as 
follows: 

Step 1 : Projection 

Projection data p % . is calculated from an estimated image I k and 
compared with the measured data p n(k+J) for all j e [0,..., M - 1] . (projection) 
15 P'j = Pn^h Vj e [0,...,M -l] 

Let A j = Kik+j)(PnVc+j) -P\- ) 
Step 2: Filtering: 
Let A 0 « A 0 

and A. = A . - J^£x*hA J=h-M-1 

20 It can be shown that one step of filtered SART gives the same result as 

M steps in ART. 

As may be taken from Step 2, the filtering is performed such that a 
product of a projection of a current angle and an accumulation of backprojections of 
preceeding angles is subtracted from the difference image A y . Advantageously, as 

25 already indicated above, this allows to filter out influences of other projections on the 
current projection. 

The above filtering step involves operations of the type P.Bj . It should 

be noted that this is a mapping from projection space into projection space. The 
combined operation can be expressed analytically and discretized in a second step. This 
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means, the backprojection and projection operations in Step 2 do not have to be carried 
out as such, in contrast to Steps 1 and 3. Only the much simpler combined operation 
P.Bj has to be carried out instead. 

Depending on the system geometry, other simplifications may apply, for 
5 example, in the case of a CT system, it is obvious that P { Bj depends only on i-j. 

Figs. 2 to 4 show images of a FORBILD thorax phantom reconstructed 
with ART, SART and filtered S ART according to the present invention (hounsfield 
units, level = 0, window = 400). All images are the result of one iteration with a 
constant value of X = 1 . Fig. 2 shows a thorax phantom reconstructed with one iteration 

1 0 of ART. Fig. 3 shows a comparison of SART (left) and filtered SART (right) for M = 
32. Fig. 4 shows a comparison of SART (left) and filtered SART (right) for M= 64. 

As may be taken in particular from the SART images on the left sides of 
Figs. 3 and 4, the image quality of the SART images is worse than the image quality of 
the ART image depicted in Fig. 2. This is due to the conservative up-date weighting of 

15 1/M On the other hand, as may be taken from the filtered SART images on the right 

side of Figs. 3 and 4, the filtered SART images have substantially the same or may even 
have an improved quality in comparison to the ART image depicted in Fig. 2. 

As indicated above, the above described image processing device and 
method of performing an iterative data reconstruction may, in particular, be applied in 

20 computed tomography. However, it may also be applied in nuclear imaging or X-ray 
imaging. 



